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A Model of the Heliosphere with Jets 

J. F. Drake^ M. Swisdak^ M. Opher^ 

ABSTRACT 

An analytic model of the heliosheath (HS) between the termination shock 
(TS) and the heliopause (HP) is developed in the limit in which the interstellar 
flow and magnetic held are neglected. The heliosphere in this limit is axisym- 
metric and the overall structure of the HS and HP are controlled by the solar 
magnetic held even in the limit in which the ratio of the plasma to magnetic held 
pressure, /3 = SttP/R^, in the HS is large. The tension of the solar magnetic 
held produces a drop in the total pressure between the TS and the HP. This 
same pressure drop accelerates the plasma how downstream of the TS into the 
North and South directions to form two collimated jets. The radii of these jets 
are controlled by the how through the TS and the acceleration of this how by 
the magnetic held - a stronger solar magnetic held boosts the velocity of the 
jets and reduces the radii of the jets and the HP. Magnetohydrodynamic (MHD) 
simulations of the global helioshere embedded in a stationary interstellar medium 
match well with the analytic model. The results suggest that mechanisms that 
reduce the HS plasma pressure downstream of the TS can enhance the jet outhow 
velocity and reduce the HP radius to values more consistent with the Voyager 1 
observations than in current global models. 

Subject headings: ISM: jets and outhows - Stars: jets - Sun: heliosphere - Sun: 
magnetic helds 


INTRODUCTION 


The historically accepted shape of the heliosphere is that of a comet-like object with a 
long tail that is dragged downstream by the flow of the local interstellar medium (LISM) past 
the sun fjParkerlll96ll: iBaranov fc Malamalll9931) . These early pictures, however, were based 
on a hydrodynamic description of the solar outflow - the solar magnetic held was assumed to 
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play a negligible role in the overall structure of the heliosphere and its interaction with the 
solar wind. Computational models based on the magnetohydrodynamic (MHD) equations 
included the solar magnetic held as we ll as that of the interstellar medium and also produced 


a heliosphere with a comet-like shape flOpher et ah 

20061 

Pogorelov et ah 

2007; 

Opher et ah 

2009; 

Washimi et ah 

2011; 

Pogorelov et ah 

20131: f 

)pher & Drakel 

2013h. 


On the other hand, the measurements of energetic neutral atoms (ENAs) by IBEX 
and CASSINI produced some surprises. These ENAs travel long distances through the 
heliosphere without being inhuenced by the ambient magnetic held and therefore yield in¬ 
formation about the large-scale structure of the heliosphere. The CASSINI ENA huxes from 
the direction of the nose and the tail were comparable, leading the CASSINI observ ers to 
conclude that the heliosphere was “tailless” flKrimigis et al.ll2009l: iDialvnas et al.ll2013l) . The 
IBEX observations from the tail revealed that the hardest spectrum of EN As were localized 


in tw o lobes at high latitude while the softest spectra were at low latitudes flMcComas et ah 


20131 1. 


Recent MHD simulations using a monopole model for the solar magnetic held, designed 
to reduce the numerical dissipation of magnetic energy that arises from a conventional dipole 
model, revealed that the solar magnetic held was str ong enough to coll imate the solar wind 
into a pair of jets that how to the North and South fjOpher et al.ll20151) . These jets bend in 
the direction of the tail, pushed by the how of the LISM. The interstellar rather than the 
solar wind plasma hows between these jets in the equatorial region downstream. Such bent 
jets h ave been seen in protost ellar systems flFendt &: Zinneckerl Il998l: ICueth fc Guilloteau 
1999) and clusters of galaxies ( Owen fc Rudnick 19761) . Astrophysical jets around 


massive 


black holes are thought to be driven by centrifugal for ces that sling the plasma along a 
rotating helical magnetic held (IBlandford fc Pavnelll982|) . However, the jets in the case of 
the helio sphere are driven in t h e region downstream of the TS as was pr oposed for the Crab 
Nebula flBegelman fc Lil Il992l: IChevalier fc Luol Il994l: iLvubarskvi 120021 1. In this region of 
subsonic how, the magnetic tension (hoop) force is strong enough to collimate and drive the 
wind. 


MHD models of the global heliosphere are complex and the mechanisms that control the 
shape of the HP, the thickness of the HS, the structure of the heliospheric jets, including the 
driver for the outhow, remain uncertain. The Voyager 1 observations have revealed that the 
thickness of the HS is around 30Af/, which is substantially thinner than expected from the 
global simulations. We present an analytic model of the heliosphere outside of a spherically 
symmetric TS where we neglect the ambient flow and magnetic held of the LISM. Taking the 
resulting heliosphere as axisymmetric and the hows within the HS as subsonic, we obtain 
the pressure and magnetic held structure of the HS along with the radius Vhp of the HP. 
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The overall shape of the HS takes the classic form of an astrophysical jet: the flows through 
the TS are accelerated to the North and South by the solar magnetic held. The heliopause 
radius is determined by continuity: the plasma how through the TS must balance the outhow 
through the jets. We present parallel global MHD simulations in the limit of zero magnetic 
held and how in the LISM which support the analytic model. 


One reason the inhuence of the solar magnetic held on the structure of the heliosphere 
is often neglected in the literature is because the pressure of the ambient plasma is large 
compared with that of the magnetic held - /3 = STmT/B"^ ~ 10 just downstream of the TS. 
We show, however, that the total plasma pressure does not control either the hows in or the 
thickness of the HS. The overall pressure in the HS is balanced by the pressure in the LISM. 
It is the tension force of the H S magnetic held that controls the pressure diherence between 
the TS and HP f AxfordI 1972 ). To the North and South there is no tension force and this 
same pressure diherence drives the axial how of the heliospheric jets. Thus, it is ultimately 
the solar magnetic held that controls the large-scale structure of the HS. 


2. Analytic Model of the Heliosheath and Heliopanse 

We consider a simple axisymmetric system in which there is no LISM how or magnetic 
held and the LISM is specihed by its ambient pressure Plism- We write down the steady- 
state MHD equations, including continuity, pressure, momentum and magnetic held, 

V • nV = 0, (1) 

V ■ = 0, (2) 

( td2 \ td2 

P + —-Vr, (3) 

Svr / dvrr 

V X (V X B) = 0, (4) 

where r is the radius in cylindrical coordinates, T is the ratio of specihc heats and B is in the 
azimuthal direction. These equations are solved in the HS with boundary conditions on the 
density, pressure, magnetic held and how given just downstream of the spherical TS which 
has a spherical coordinate radius F^. At the TS we assume that the azimuth how is zero 
and since there are no forces in the (j) direction (Eq. ([3])), we can take Vj* = 0 everywhere. 
Thus, from Eq. ([T]) we can write 


nV = V(/> X Vi/j, 


(5) 
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where "0 is the stream function for the particle flux. Taking B = rBVcj) with V from Eq. (1^ . 
we can reform Faraday’s law in Eq. (|1]) as 


X V 



0 . 


( 6 ) 


This equation has components only in the 0 direction so taking the dot product of this 
equation with V0 yields the constraint that B/nr is constant along streamlines or 


B/nr = fi'ip), 


(7) 


where the function / is only a function of ^|J. The form of / can be determined by the 
boundary conditions along the TS. Downstream of the TS we take the flow to be normal to 
the shock, in the radial (in spherical coordinates) direction with a constant value Vg with 
the density also constant. So, from Eq. (jS]) we find dip/89 = —R’^ngYsSmO or 

ip = ngVgRlcosOg, (8) 


where 9g is the polar angle in spherical coordinates at the shock. Thus, the variation of 
ip along the shock is known. Similarly, we know from the solutions of the Parker spiral 
magnetic held that along the shock B = Bgsmdg so B/nr = Bg/ugRg is a constant and so 
is / = Bg/ugRg. Throughout the HS we have 


B _ Bg 
nr ngRg 


(9) 


Turning to the pressure equation and using Eq. (|5]), we hnd V0 x Vip ■ V{P/{nF) = 0 so 
P/n^ is also constant along a streamline and therefore a function only of ip. As before, we 
can evaluate it along the TS where it is given by Po/n/.. Thus, 


P = Ps 




( 10 ) 


Thus, both P and n in the HS are linked to B and r. 


We now focus on the high f] limit, which is most relevant to the HS, where the pressure 
associated with interstellar pickup ions dominates the magnetic pressure. Specihcally, we 
take to be a large parameter. In addition, since the flows are subsonic downstream of the 
TS, inertial forces are also small. We can therefore write the plasma pressure in a series 
Po + Pi + ... with Pi ~ B'^/Sti ~ MnV^/2 -C Pq. Thus, to lowest order Eq. ([3]) becomes 


0 = -VPo(P,r), 


(11) 




5 


where Pq is an explicit function of B and r through Eq. flTOl) . Equation ffTOj) requires that 
to lowest order the pressure in the HS is constant everywhere and is given the value Pg at 
the TS. The density is also a constant, rig. Since P is linked to B and r through Eq. flTUD 
the constancy of Pq requires that 


B 


Bg — 
R. 


( 12 ) 


SO that B increases with radius outside of the TS f)Axfordlll972l: IChevalier fc Luolll994l) . At 
hrst order we include the inertial terms and magnetic field in the momentum equation, which 
becomes 

P?r2 \ PV 


MugVo ■ Wo = -V Pi + 


871 R^g) 




Vr, 


(13) 


where n has been replaced by Ug and Vq = V0 x V'lp/ng. 


Before discussing the flows in the HS, we consider the weak flow limit of Eq. flT^ so that 
the inertial forces in the radial direction can be discarded. In this limit the magnetic tension 
force in Eq. flT^ causes the plasma pressure and total pressure to decrease with radius. This 
limit is artihcial for the HS since the flow Vg downstream of the TS is comparable to the 
Alfven speed and the associated radial inertial forces are comparable to the magnetic forces. 
Nevertheless, this limit illustrates how the pressure in the HS varies. We do not require zero 
pressure gradient in ; 2 . The pressure drop from the TS to the LISM is balanced by magnetic 
tension in the radial direction but the same pressure drop also develops from the equator to 
the outflow jets to the North and South. This pressure gradient along z drives the outward 
flows associated with the jets. Thus, we integrate Eq. flT^ from the TS outwards to obtain 
an explicit expression for Pi, 


Pi(r-) 


——ir 

47rP2 ^ 


Rl sin^ Og), 


(14) 


where 6s{z) is the value of 6 at the TS and is dependent on z. Pressure balance across the 
HP, which requires that P{rhp) + B‘^{rhp)/8n = Plism, then yields an explicit expression for 


^hpj 


' hp 


SttAP 


+ 2sin^ 6g, 


(15) 


where AP = Pg — Plism- At this stage in the calculation the pressure difference AP remains 
undetermined. We will show, however, that the requirement that the mass flow into the HS 
across the TS balance that out the two jets constrains the pressure difference. In Fig. [T]we 
show 2-D plots of the plasma pressure, the magnetic pressure and the total pressure in the 
r, 2 ) plane in the HS. The data is shown for SttAP/B^ = 2, which as shown later is the upper 
limit on AP. The inner boundary of the data shown is the TS and the outer boundary is 
the HP. The radius of the HP peaks at the midplane and falls off with z until z > Rg, where 
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it remains constant, forming the Northward jet. The plasma and total pressure decrease 
with radius r while the magnetic pressure increases with r. The total pressure falls off with 
distance from the equator until it approaches a constant value in the jet. We will show that 
this pressure difference, which is a consequence of magnetic tension, drives the jet outflow. 
Along the axis P remains constant at Pg. This shape was obtained by neglecting the radial 
plasma inertia, an assumption which breaks down where the straight portion of the HP in 
Fig. [1] intersects the curved portion at z/Rg = 1. The sharp kink in the HP is not real and is 
not seen in the MHD simulations discussed later. In Fig. HJcuts along r of the total, plasma 
and magnetic pressures from the data of Fig. [T]are shown at the equator in (a) and across 
the jet in (b). 

We now discuss the flows driven in the HS. We derive a Bernoulli-like equation by taking 
the dot product of Eq. ffT5]l with Vq and integrating along the streamline. 


2.^2 


1 


= -Mugv: + 


Bl sin^ Og 
An 


(16) 


Unfortunately, this equation does not enable us to calculate Vq throughout the HS because 
it requires that we know the trajectory of a stream line within the HS to link the local radius 
r with its position at the TS where 6g is known. Along the axis where = 0 Eq. 0161) gives 
kb = Vg, since the pressure is constant and B is zero. Similarly, the velocity at the jet radius 
r^ei, can be calculated since at that location Og = n/2. We hnd V^iTjet) = where 

= Bl/AnMug. Thus, the increase in the jet velocity above Vg is linked to the Alfven 
velocity based on the magnetic held strength Bg at the TS. More generally, we can calculate 
Vb(r) across the jet radius by noting that within the jet 


Uo = Uo, = 


1 dijj 
UgT dr 


(17) 


From Eq. ([8]) sin 6s can be written in terms of "0 so we are left with a single equation for ip 
across the jet. 


1 





U/ + 2c 


2 

As 



r \ 

n^sRtVgy ’ 


(18) 


where ip varies from ngVgR^ at the jet axis to 0 at the HP. The equality of the particle huxes 
through the TS and jets requires that Eq. flTSD produce the requisite jump in ip across the 
jet. Equation (1T8|) can be simplihed by dehning an angle variable cos 6* = ip/ngVgnl, 




89 




(19) 


where 9 varies from 0 at the jet axis to n 12 at the HP. This equation can be integrated 
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directly to obtain the jet radius rjet, 



^ 2 tan ^{\/2cas/Vs) 
^ V2cas/V, 


( 20 ) 


The jet radius is a maximum for cas ^ K when the jet outflow velocity is given by Vg. In 
this limit the conservation of particle flux reduces to the jet cross-sectional area being equal 
to the TS area or rjet = V^Rs- With increasing cas the outflow velocity of the jet increases 
and Tjet decreases. For cas S> Vg, rjet oc Rg^JVs/ cas- An expression for the pressure jump 
AP between the TS and the LISM can be calculated from Eq. flT^ . which is exact in the 
jet where Vq,. = 0, 

AP = , (21) 

Stt ^CAs/Vg 

The pressure jump is a maximum when cas is small and decreases with increasing cas- The 
dependence of r^et and AP are shown as functions of CAg/Vg in Fig. [31 From Eq. flTS]) 
therefore also decreases with increasing CAs/Vg. 


3. Global MHD Simulations 

We have carried out MHD simulations of the global heliosphere without an interstellar 
wind and magnetic held. Our model is based on the 3D multi-huid MHD code BATS-R- 
US. It envolves one ionized and four neutral H species as well as the magnetic held of the 
sun. We used a monopole conhguration for the solar magnetic held to eliminate artihcial 
reconnection across the heliospheric current sheet. The basic parameters of the simulation 
are the same as those described in Opher et al. 2015. The computational grid was ±3000 AU 
in each direction. Parameters of the solar wind at the inner boundary at 30 AU were: vgw = 
417fcm/s, nsw = 8-74 x 10“^cm“^, Tsw = 1-087 x 10^77 and the Parker spiral magnetic 
held with a radial component Bsw = 7.17 x 10“^nT at the equator (with an azimuthal 
component = 0.22nT). The solar wind how at the inner boundary is assumed to be 
spherically symmetric and the magnetic axis is aligned with the solar rotation axis. For the 
LISM we assume Tism = 6519 K while the plasma density was raised to 0.483/cm^ to make 
up for the absence of pressure associated with the interstellar magnetic held. The number 
density of H atoms in the interstellar medium is uh = 0.18 cm~^ and the temperature is the 
same as for the interstellar plasma. The z-axis is parallel to the solar rotation axis. The grid 
has cells ranging from 2.93 AU at the inner boundary to 187.5 AU at the outer boundary. The 
simulation had a resolution of 3.0HP between z = ±750HP and x = ±305HP; y = ±400HP, 
encompassing the entire HS. The run was stepped forward for 3061 years. 





In Fig. m we show in yellow the surface of the HP as dehned by InT = 13.9. The 
simulation reveals jets to the North and South as in the analytic model. The HP bulges at 
the equator as in the model. The gray lines are the solar magnetic held. Shown in Fig. |5]in 
the x — z plane are the plasma pressure in (a), the magnetic pressure in (b) and the speed and 
streamlines in (c). As in the model, the plasma pressure decreases with cylindrical radius 
r away from the jet axis while the magnetic pressure increases with r and the strongest 
magnetic helds are in the equatorial region just upstream of the HP. The streamlines reveal 
the North and South directed outhows that make up the jets. The HP boundary does not 
reveal the sharp indentation seen in the model. Finally, in Fig. |2]J^c) we show cuts of the 
total pressure (solid), the plasma pressure (dotted), the magnetic pressure (dashed) and the 
magnetic held (dot-dashed) in cuts along r at the equator from just upstream of the TS 
to past the HP. The pressures have been normalized to B to Bg and r to Rg with 

Rg = 135AH taken to be the location of the maximum of Ppiasma- The cuts are in remarkable 
agreement with the cuts from the model in (a). In the simulation cas/Vs = 1-1, which from 
Fig. [3], yields SnAP/Bg = 1.3 compared with the measured value of 1.5 from Fig. [2](c). For 
StiAP/B] = 1.5 Eq. (IT^ yields rhp/Rs = 1-9 at the equator, essentially identical to the 
HP radius in the cuts, and the ratio of the HP radius at the equator to that in the jet is 
0.64 (Eq. 0151) I compared with the reasured value of 0.66. Finally the measured particle flux 
through the TS is the same as that out the jets. 


4. Discussion and Conclusions 


We have explored the structure of the HS and HP when the interstellar flow and magnetic 
held are neglected and the system can be treated as axisymmetric. We show that even in the 
limit in which P ^ in the HS the magnetic held controls the large-scale structure of 

the HS and drives Northward and Southward directed jets. To lowest order the pressure in 
the HS is balanced by the pressure in the interstellar medium. The magnetic held controls 
the pressure variation within the HS and re-directs and boosts the how across the TS to the 
North and South to form heliospheric jets. The radial distance from the TS to the HP and 
the jet radii are controlled by the requirem ent that t he plasma howing into the HS across 
the TS hows outwards in the jets (see also lYul (119741) 1. For very weak magnetic helds the 
jet outhow velocity is the same as the velocity Vg downstream of the TS. In this limit the 
total cross-sectional area of the jets is equal to the area of the TS and rj^t = \/2Rg. With 
increasing magnetic held strength the jet outhow velocity increases and the radii of the HP 
and the outhow jet decrease (Eq. fl20l) and Fig. [3](a). 


The global MHD models of the heliosphere flMalama et al.ll2006l: iPogorelov et al.l 120071: 












9 


Opher et al.ll2009l) produce HS thicknesses that are around 50 — 70^4f/, substa ntially larger 


than the value of SOAf/ determined from Voyager I’s crossing of the HP in 2012 fjStone et ah 


2013|). The results here suggest that mechanisms that increase the jet outflows will reduce the 
HP radius. Pressure reductions in the downstream region associated with thermal conduction 
or other mechanisms might produce such enhanced flows. 

Th ere is evidence fro m both the present MHD simulations and those carried out the 
earlier (jOpher et al.ll2015h that the jets are subject to l arge-scale instabilities. T he resulting 
turbulence might be a driver of anomalous cosmic rays flStone et al.ll2005l. l2008l) . High time- 
resolution ENA measurements might be able to establish the existence of the heliospheric 
jets and associated turbulence. For a jet radius of around 140A7/ and an Alfven velocity of 
around lOO/cm/s, the Alfven transit time is around 10 years. The jet turbulence might, of 
course, cascade to smaller scales so the relevant time scales could be shorter. 


This work has been supported by NASA Grand Challenge NNX14AIB0G and NASA 
awards NNX14AF42G, NNX13AE04G and NNX13AE04G. The MHD simulations were car¬ 
ried out on Pleades at the NASA Ames Research Center under the award SMD-14-4986. 
We acknowledge support from the International Space Science Institute for the team “Fac¬ 
ing the Most Pressing Challenges to Our Understanding of the Heliosheath and its Outer 
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Fig. 1.— 2-D images in the HS of the plasma pressure in (a), the magnetic pressure in (b) 
and the total pressure in (c), all normalized to Pbs = where Bg is the magnetic held 

just downstream of the TS at the equator. The images are for AP = Pg — Plism = 2i?^/87r 
and /9s = Ps/Pbs = 12, where Pg is the plasma pressure downstream of the TS and Plism 
is the pressure of the LISM. 
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Fig. 2.— Cuts through the data of Fig. [T] of the total pressure (solid), plasma pressure 
(dotted) and magnetic pressure (dashed) versus r at the equator in (a) and within the jet in 
(b). In (c) cuts from the MHD simulation of Fig. |4] along the equator. 















LISM 



"ig. 3.— The jet radius rj^t and the pressure jump between the TS ai 
— Plism^ versus the Alfven speed downstream of the TS at the e( 
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Fig. 4.— The heliopause as defined by InT = 13.9 from an MHD simulation embedded in 
an ambient interstellar medium with no mean flow and zero magnetic field. The gray lines 
are the solar magnetic field with the TS visible as a disc. 
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Fig. 5.— The plasma and magnetic pressures (Pa) in (a) and (b), and the plasma speed 
{km/s) and streamlines in (c). All in the x — z plane through the center of the heliosphere 
from the simulation in Fig. HJ 


















